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Abstract 
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integral's domain of integration to the first Gribov region when the Landau gauge is considered. 
We identify a local polynomial action and prove the renormalizability to all orders of perturbation 
theory by employing the algebraic renormalization formalism. Furthermore, we also pay attention 
to the breaking of the BRST invariance, and to the consequences that this has for the Slavnov- 
Taylor identity. 
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1 Introduction 



It is a well known fact that SU(N) Yang-Mills gauge theories, described by the EucUdean action 

SYM = \j d\F^F^^, (1.1) 

with the gauge potential and 

f;, = a^e - ^vAi + gf'A'^A'^ , (1.2) 

the field strength, whereby D|^^ is the covariant derivative in the adjoint representation, given by 

< = aX'-^r'"4> (1.3) 

are asymptotically free at very high energies |[Il|2l[3l|4l. The coupling constant is sufficiently small 
to allow for a perturbative description, with asymptotic degrees of freedom given by massless gauge 
bosons. We shall not consider fermion matter in this paper, however the same conclusion holds for 
quantum chromodynamics (QCD), written in terms of gluons and quarks. When we pass to lower 
energies, the coupling constant begins to grow, and perturbation theory starts to loose its valid- 
ity. At still lower energies, the situation becomes dramatically different, as perturbation theory now 
completely fails, and confinement sets in, meaning that the elementary field excitations are no longer 
physical observables, but become confined into colorless states. The hadrons constitute the physical 
states of QCD. 

A satisfactory understanding of the behaviour of Yang-Mills theories in the low energy regime is yet 
unavailable. Due to the large coupling constant, nonperturbative effects have to be taken into account. 
The introduction of condensates, which are the (integrated) vacuum expectation value of certain local 
operators, allows one to parametrize certain nonperturbative effects arising from the infrared sector 
of e.g. the theory described by (11.11 ). Condensates give rise to power corrections, a phenomenon that 
can be handled using the Operator Product Expansion. Clearly, these power corrections correspond to 
nonperturbative information in addition to the perturbatively calculable results. If one wants to con- 
sider the possible effects of condensates on physical quantities in a gauge theory, only gauge invariant 
operators are relevant. The most famous example is the dimension 4 gluon condensate (^a^F^^'j, giving 
rise to ^ power corrections in QCD. The SVZ (Shifman-Vainshtein-Zakharov) sum rules ||5] can then 
be used to relate the condensates to observables, and hence one may obtain certain phenomenological 
estimates for e.g. (ttv^^/v)- This approach allows for a study of at least some aspects of QCD in an 
energy regime in between the confined and perturbative zone. One can still use perturbation theory 
using the quarks and gluons as effective degrees of freedom due to the lack of explicit knowledge of 
the correct physical degrees of freedom, but the results get adapted by nonperturbatively generated 
condensates. 

In this paper, we shall introduce an action that can serve as a starting point for investigating some 
nonperturbative effects in gauge theories. These nonperturbative effects arise from 2 premises: the 
Gribov problem in fixing the gauge freedom, and the possibility of a dynamical mass generation. 

First of all, we shall be concerned about fixing the gauge. If we want to perform any kind of calcu- 
lations, we must reduce the enormous gauge freedom of (11.11 ). encoded in the local gauge symmetry 
generated by 

doA'l, = -of to'' , with (0^ arbitrary , (1.4) 
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to a global one by a suitable gauge fixing condition, say ^ (A^) = 0. In principle, by imposing a 
gauge condition one should select a single representative A* from the gauge orbit A^, where U is a. 
generic SU (N) gauge transformation. Unfortunately, it was shown that is impossible to uniquely fix 
the gauge, a problem related to the complicated topology of the space of gauge orbits |6|. 

Seminal work on the existence of gauge copies was done 3 decades ago by Gribov in |7 |. This paper is 
not the place to give a complete overview of the ambiguities arising when a gauge fixing is performed, 
we therefore kindly refer to Gribov 's original paper or to the available literature, such as |8J, which 
contains many examples and references. In particular, in [7|, it was pointed out that, in the Landau 
and Coulomb gauges, the existence of zero modes in the Faddeev-Popov operator gives rise to gauge 
copies. Using the Landau gauge, 

d,A, = 0, (1.5) 

one finds that a gauge equivalent configuration Aj[,, connected to A^, via ( 11.41 ). also obeys 3pAj[, = 0, 
when 

fM-"V = 0, (1.6) 
where 5W "'^ denotes the Faddeev-Popov operator 

^"' = -d,{d,?>'"'-gr''Ai,). (1.7) 

The existence of the Gribov copies implies that the domain of integration in the path integral has to be 
further restricted in a suitable way. Following Gribov, it seems logical to restrict to the region Q. with 
corresponding boundary dQ., which is the first Gribov horizon, where the first vanishing eigenvalue 
of the Faddeev-Popov operator (11.71 ) appears Q. Within the region Q. the Faddeev-Popov operator 
is positive definite, i.e. M"'^ > 0. Quite obviously, this restriction to the first Gribov region can be 
motivated only if every gauge orbit passes through it. It was shown by Gribov that this is certainly the 
case for gauge potentials "sufficiently close" to the boundary ^^l 171, whereas the proof for general 
configurations was presented in @|. Nevertheless, we should also mention that the Gribov region 
itself is also not free from gauge copies f9j UHl [HI [El- To avoid the presence of these additional 
copies, a further restriction to a smaller region A, known as the fundamental modular region, should 
be implemented. Nevertheless, the implementation of the restriction of the domain of integration to A 
proves to be a quite difficult task which, to our knowledge, has not yet been accomplished. Recently, 
it has been argued that the additional copies existing inside D. might be irrelevant when computing 
expectation values, meaning that averages calculated over A or H should give the same value [.131 . 

Using a semiclassical argument fj), Gribov implemented the restriction to the region Q.. Essentially, 
his argument relied on the fact that the (Fourier transform of the) inverse of the Faddeev-Popov oper- 
ator, which is nothing else than the ghost propagator, encounters no poles elsewhere than at the origin 
= 0. This amounts to say that the operator !M "'^ itself does not vanish, except at the horizon. By 
using this "no pole condition", we are assured that the considered gauge potentials remain inside the 
first Gribov regioifl, and as a such at least the set of copies obtained via (11.6b is already excluded from 
the game. 

This restriction has many important consequences for the infrared behaviour of the propagators. The 
gluon propagator turns out to be suppressed in the infrared, while the ghost propagator gets enhanced 
fTl . Moreover, it can also be shown that the gluon propagator exhibits a violation of positivity in 
its spectral density representation, a sign that the gluon cannot be a physical observable anymore, 

'Albeit that this restriction cannot be implemented exactly, but only in an order by order expansive way. 
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see lfT4l [T5l [T6l [TTl and references therein. It is interesting to mention that lattice simulations of the 
Landau gauge propagators have revealed evidence for this suppression, respectively enhancement, 
see e.g. |[T8l [111 |2ll l22l |23l IM ISll |26|. Another consequence of the Gribov restriction is the 
"infrared freezing" of the strong couphng constant, i.e. Ois{p^) tends to a constant as goes to zero, 
see lfT4l |27]| and references therein. Again, this behaviour is in qualitative agreement with lattice data 
II22I l23l l24l as well as with the results obtained from the analysis of the Schwinger-Dyson equations 
llll |29l [301 [3B [321 [31 [M [351 [Ml ■ 

It might be clear that the restriction to the Gribov region Q. could be of great relevance for a better 
understanding of the infrared region of gauge theories. This belief is further supported by the Kugo- 
Ojima confinement criterion |[37l which, in the case of the Landau gauge, turns out to rely on a ghost 
propagator diverging stronger than [38]. This feature is also present when the restriction to the 



p- 



Gribov region is implemented, yielding in fact a ghost propagator developing a singularity. 

An important progress on the restriction to the Gribov region Q. was accomplished by Zwanziger in the 
papers [ 39^,'40'|. The restriction to Q. was implemented through the introduction of a nonlocal horizon 
function appearing in the Boltzmann weight defining the Euclidean Yang-Mills measure. According 
to |[39l[40l . the starting Yang-Mills measure in the Landau gauge is given by 

dfjy = DADcDcDbe ('^+^^) , (1.8) 

where the starting action is 

5 = 5ym + %, (1.9) 

with Sgf the gauge-fixing action given by 

% = I d\ (b^'d.Al + 7^ d.Dfc'') , (1.10) 

where the auxiliary field b"^' is a Lagrange multiplier enforcing the Landau gauge (11.51 ). (c'^,c'^) are the 
Faddeev-Popov ghost fields, and 

H = I d'xh{x) =g^l d'xr'^A'^i^-Yf'^, (1.11) 

is the so-called horizon function, which implements the restriction to the Gribov region Q.. We recog- 
nize that H is nonlocal. The massive Gribov parameter y is fixed by the horizon condition 

{h{x))=4{N^-l) , (1.12) 

where the expectation value {h{x)) has to be evaluated with the measure (11.81 ). To the first order, the 
horizon condition (11.121 ) becomes, in d dimensions, 

_N{d-\) 2 [ d^'q \ 
^~ 4 ^ J {2n)W + 2NgY ■ 

This equation coincides with the original gap equation derived by Gribov for the parameter y j'T^. 

We shall rely on the path integral formalism, so that we can localize the horizon fuction (11.111 ) by 
means of a pair of complex bosonic vector fields BOl . according to 



(1.14) 
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where the determinant, (det5Vf takes into account the Jacobian arising from the integration over 

((^^^^F;'),and 

f = D{N^-\)=A{N^-\), (1.15) 

with D = 4 the dimension of the Euclidean space time, and N the dimension of the gauge group. 
This determinant can also be localized by means of suitable anticommuting complex vector fields 
((o;;^ namely 



(1.16) 



(det mY = j DToD(o exp - J A (^-W"^ "^(O^ 
Henceforth, the nonlocal action Su is transformed into a local one given by 



where 
and 



Sy ■ 



J A [<^1"m - , ( i . i s) 

jd\gf''{r;-f;)Ai. (1.19) 

As it was shown in |[T4l [39l l40l |4TI . the resulting local action turns out to be renormalizable to all 
orders of perturbation theory. This is a point of great importance, as it allows for a consistent and 
order by order improvable framework to calculate relevant quantities when the restriction to the Gribov 
region Q. is taken into account. 

A second point that motivated this paper is the issue of the dynamical mass generation in gauge 
theories, and related to it that of the ^ power corrections. A few years ago, in a series of papers, 

Zakharov et al. questioned the common wisdom that ^ power corrections cannot enter gauge invariant 
observables, as local gauge invariant operators of mass dimension two do not exist. This is a reflection 
of the fact that one cannot add a renormalizable mass operator for the gauge fields to the Yang-Mills 
action, at least not when the Higgs mechanism and associated symmetry breaking are not considered. 
However, by using QCD sum rules, it was advocated in [42] that an effective gluon mass could account 
for the \ corrections, leading to an acceptable phenomenology. The underlying condensate was 
proposeJ^ to be the gauge invariant quantity ||43l |S1 



which originates from a highly nonlocal operator, since 1431 and references therein] 

+ 0{A^) 



a2, _ 1 / ^4 



and therefore it falls beyond the OPE applicability. The interest was especially focused on the Landau 
gauge, since then the operator A^j^^ reduces to the local quantity A^. An effective potential for (A^) was 
calculated up to two loops in H6ll47]| . giving evidence for a nonvanishing condensate and consequent 
effective gluon mass oc (A^). Determining a sensible effective potential for a local composite 
operator (LCO) is a nontrivial task, but nevertheless it was dealt with in p6l based on the method 
developed in BSll . The renormalizability of the so-called LCO method was proven in B9l to all 
orders of perturbation theory in the case of A^ in the Landau gauge. 
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Unfortunately, the gauge invariance of A^^^ is, strictly speaking, only ensured when the absolute 
minimum of A^^^ along the gauge orbit has been reached, a highly difficult task due to the presence 
of Gribov copies. Moreover, it is unclear what can be done with this operator beyond the Landau 
gauge. In other gauges, other renormalizable and condensing dimension two operators exist, but 
these are explicite gauge parameter or even ghost dependent, see ||50ll for an overview. In e.g. the 
maximal Abelian gauge, an effective mass was found for the off-diagonal gluons only f5T|, which is 
qualitatively consistent with the available lattice data Il52ll53l . 

Let us also mention that effective gluon masses have been studied in the past from theoretical, phe- 
nomenological and numerical viewpoint, see 1201 |46l |5l |55l 111 |57l |58l |59l |60l ^ from 
exhaustive list. 

Taking all this into account, it seems to be a worthy task to look for other potential candidates which 
could be at the origin of the dynamical mass generation and related ^ power corrections. It would 
also be favourable to start from a gauge invariant operator. The candidate we already investigated in 
l62l is the nonlocal operator 



O 



ah 



(1.21) 



which can be coupled to the Yang-Mills action via a nonlocal mass term 



So = -— I d\F^^ 



ah 



(1.22) 



where m is a mass parameter, and [[D ) ] is the inverse of the covariant Laplacian 



■abc~ 



(1.23) 



Analogously to what has been done in the case of the nonlocal horizon function (11.1 lb . the gauge in- 
variant mass operator ( 11.221 ) can be localized with the help of a pair of complex bosonic antisymmetric 
tensor fields in the adjoint representation P31 . {B'^^,B"i^), 



-So 



DBDB {det D^y exp<{ / A 



(1.24) 



Here we have 



D{D-\] 



(1.25) 



and, like in the case of the horizon function, the determinant, (detD^)^', can be localized using a pair 
of anticommuting antisymmetric complex tensor fields (G^v5^"iv)' according to 



(detD- 



DGDG exp 



-- I d\(Gr^,Dl^D^,'G% 



(1.26) 



Then, the action So gets replaced by its local version given by 



S Local c I o 
O — ^BG + , 



where. 



(1.27) 
(1.28) 
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and 

5,„ = ^|A(B^,-S;Jf;,. (1.29) 

We underline the fact that an initially nonlocal operator can be cast into a local form ll63l . In the case 
of A^jjj, this would not be possible, as it is a infinite series of different nonlocal operators. Once we 
arrive at a local action, we can investigate e.g. the renormalizability to all orders by means of algebraic 
methods, the canonical quantization, the explicit calculation of the renormalization factors, etc. 

The goal of this paper is to study the massive action ( 11.221 ) when the restriction to the Gribov region 
n is implemented a la Zwanziger. Since the extended action Sym + So is gauge invariant, we might 
expect that the procedure of further restricting the domain of integration will have no influence on 
the renormalizability. This will be explicitly confirmed. In a future stage of research, one can start 
searching for the value of the Gribov parameter y as well as the dynamically generated mass m. 
We remind here that the Gribov-Zwanziger action itself can also be used to mimic 4- corrections, 
as explicitly discussed in ll27l . As a future endeavour, it would be worthwhile to study physical 
correlators with our action, and find whether the Gribov and/or mass parameter m are a potential 
source of such power corrections. 

Let us return to the content of this paper, which is organized as follows. In section 2, we introduce 
all the necessary sources in order to find a suitable starting action. The set of Ward identities defining 
this action is presented in section 3, while in section 4 we compute several useful (anti-)commutation 
relations between the linearized symmetry operators. These are used in section 5 in order to construct 
the most general allowed invariant counterterm. In section 6, we confirm the renormalizability since 
we shall be able to reabsorb all the allowed counterterms in the starting action by introducing suitable 
bare quantities. In section 7 we discuss a few properties of the physical action, which is obtained 
from the starting action by setting the sources equal to their physically relevant values. The process 
of giving specific values to the sources breaks the BRST invariance. In section 8, we discuss the 
associated breaking of the Slavnov-Taylor identity, and we comment on the fact that in most cases 
this breaking becomes harmless for the identities derivable between a large class of Green functions. 
Finally, section 9 is devoted to the conclusions. 



2 Identification of the complete classical action L 

We shall start with the following local action, as it was obtained in the introduction 

e c I c I cLocal i cLocal i \ 

•JLocal — JYM + "Jgf + +"^0 • y'^-'^) 

As we wish to discuss the renormalizability, we should try to establish as many symmetries as possible. 
These symmetries can then be translated into Ward identities. As we are dealing with a gauge theory 
which is to be gauge fixed, we expect to find a BRST invariance and consequent Slavnov-Taylor 
identity. All these identities are a powerful tool in constructing the most general allowed counterterm 
ll64l . If this counterterm can be reabsorbed in the original action through the introduction of bare 
quantities, we are able to conclude that the starting action is renormalizable. If not, we could still try 
to identify a more general starting action that is renormalizable. This has been discussed in extenso 
already in fl31l62]| when analyzing the nonlocal mass term (11.221 ). 
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2.1 BRST invariance 



In order to find the BRST invariance of the resulting local theory, given by ( I2.1I ). we proceed as in 
BOl i45J and consider at first the particular case when y = m = 0, i.e., 



y=m=0 
Local 



c I c I cLocal,y=0 „Local,m=0 



— '^YM + Sgi + S(^a + Sbg ■ 

In this case, we have actually introduced nothing more than two unity factors, written as 
1 = J D^D(|)DroD(o exp 
1 = / DBDBDGDG exp 



(2.2) 



(2.3) 



Nevertheless, the action (12.21) may be written in a BRST invariant fashion. To see this, let us first 
introduce the following nilpotent BRST transformation 







-Dfc" , 


sc" 










gf^'c%, + G\ 


sG" 










gf''c'G%+Bl 








sc^ 




b\ 


sb" 




0, 


< 




<, 






0, 






-r-ab 
% ' 


-rab 




0, 






0. 



/JV 



Now, let So be the action defined by 



So = Sym + s I d\ {c"d^l + %'-^ + G;,D"jD'i,'B^ 



r\acr\cbT>b \ 
/jv) 



(2.4) 



(2.5) 



which satifies 

sSo = 0. (2.6) 

Applying the BRST transformations (12.41 ) and recalling that the Faddeev-Popov operator, M"^, is 
given by ( 11.71 ). we obtain 



(2.7) 
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Following BOl one can show that and Sj^J^^^ ^ are equivalent. More precisely, one may transform 
^Lo7^^ into Sq by performing the following shift in the variable co^'', 



CO" 



CO" 



(2.8) 



and keeping in mind that the corresponding Jacobian turns out to be field independent. Thus, the 
following equivalence holds. 



D<t>e 



-So 



(2.9) 



where <I> is a shorthand for all the fields. Now, let us reintroduce the term Sy, given by ( 11.191 ). while 
Sm remains absent. It is easy to show, using the BRST transformations (I2.4I ). that Sy may be rewritten 
as 

S, = ^ I d'x[gf''r;A'-s{gf''^'A'^)+gf^^^^^^^^ . (2.10) 

The last term can be eliminated by means of a change of variables 



co^;'' ^ CO*'' + 7^ (fW- -^f''gf''D';c" 



(2.11) 



Furthermore, we notice that, thanks to fact that the integral of a total derivative vanishes, the following 
expression for Sy holds 



.ahj^ba ^ 



Nevertheless, the action. 



So + S. 



(2.12) 



(2.13) 



is not yet BRST invariant. This point can be dealt with by means of the introduction of a pair of BRST 
doublets of local external sources BOl . {M^y,N^y) and {M"f_!^,N"!^), which transform as 



sm: 



ah 



:ab 



-n: 



ab 



sM^I^ = . 



(2.14) 



As pointed out in BOl . the introduction of these external sources allows us to promote expression 
( 12.131) to a BRST invariant action. In fact, let ^sources be the action 



which obviously satisfies 



When the sources {,M^^,m"J^,N^^,N"J^) attain their physical values EOl, defined by 



(2.15) 
(2.16) 



-m: 



ab 



n: 



■ab 



phys 
phys 



phys 



phys 



0, 



(2.17) 
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it immediately follows that 



phys 



Sy = -'fld''x[D"X"-s{Dfm';^ 



(2.18) 



One sees thus that the use of the external sources (M^v'^i^vi^^^w'^v) enables us to introduce an 
extended action Zq, given by 

'-'sources ; 

(2.19) 

which enjoys the important property of being BRST invariant, 

ilo = , (2.20) 

while reducing to expression ( 12.131 ) when the sources attain their physical values, given by ( 12.171 ). 
Recapitulating, we have rewritten 

J DADcDcDbe-^+^" = J DADcDcDbD^Di^(oD(oDe^° . (2.21) 

It is then easily shown, upon combination of ( |1.11| ). (I1.12| ). ( 12.131 ) and ( 12.211 ) that the horizon condition 
is implemented by requiring that 



9y2 



0, withy^^o, (2.22) 



whereby Fgz is the Gribov-Zwanziger effective action defined by 

= J DADcDcDbD^Di^DToDoje^^ . (2.23) 
To continue, let us analyze the term S,,,, given by (fL29l ). This term is, just as S^f""^^ in ([127]) . left 



li45J 






= -DfQ\ 



















(2.24) 

where is the parameter of the gauge transformation, but it is not invariant by the BRST transfor- 
mations (12.41) . This problem can be solved in a way equivalent as done in the case of Sy. This time 
we will introduce a pair of BRST doublets of external sources, (t/ap^vi^ap^iv) ^rid (t/ap^ivi^ap^iv)' 
transforming as 

sUap^v = Voip^v , ■5VaPA(V=0. (2.25) 
Hence, by considering the following term 

•^sources — ^ J ^ i^a^^/LivGa^ — ?7ap;uv^ap)^^v ) 
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the term S'^^^^^^^ reduces to Sm of ( I1.29I ) when the sources (t/ap^vi ^^aP/jV) Vappyj^ap^iv) attain the sub- 
sequent physical values 



^ap;uv 



phys 



phys 



im 

phys 2 

= 0. 

phys 



(5a^/5pv - 6av5p^) 



These sources enable us to define an action Zi as 

II = lo + 5'sou 

in such way that 

sli = . 



(2.26) 

(2.27) 
(2.28) 



2.2 The global U (/) and U (/) symmetries 

In addition to the BRST invariance the action Ei displays global symmetries U (/), / = 4{N^ — 1) and 
U if), f = 6, respectively expressed by 



/ a™ 



-ch SZl ^„ §Ei 



5m: 



and 



51, 



u 8y,i \ 

= 0, 



SZi SZi 



5£i — c-fo 5£i 



(2.29) 



5£i 



ou, 



aP 



^JVCp 



-u 



51, 



/ivop"; 



51, 



apcp 



6£, 



OV'ccpap 



0, 



(2.30) 



O^^apop w^ivcp 

The presence of the global invariances U (/) and U (/') means that one can make use B0ll45l of the 
composite indices / = (a,/^), 7 = 1,...,/, and / = (/i,v), / = 1, . . . ,/'. Specifically, setting 



and 



fM""^ M"'' N"'' N"'') 



— (R" r" 



(2.31) 



(2.32) 



we rewrite Ei as 



(2.33) 
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A 


b 


c 


c 






CO 


CO 


B 


B 


G 


G 


dimension 


1 


2 


2 





1 


1 


1 


1 


1 


1 


1 


1 


ghost number 








-1 


1 








1 


-1 








1 


-1 


2/-charge 














1 


-1 


1 


-1 














2j/ -charge 


























1 


-1 


1 


-1 



Table 1 : Quantum numbers of the fields 





MMNNVVUU 


dimension 
ghost number 
2/-charge 


2 2 2 2 1 1 1 1 

1 -1 1 -1 

1 -1 1 -1 
1 -1 1 -1 



Table 2: Quantum numbers of the sources 



For the symmetry generators, we have 



and 



C ^ ? 

M'\^^+N",——-NIj^^ ] , (2.34) 



By means of the trace of these operators the /(/)-vaIued fields turn out to possess an additional quan- 
tum number, displayed in Tables [T] and |2j together with the dimension and the ghost number. 



2.3 The complete classical action E 

We proceed by establishing the complete set of Ward identities which will enable us to analyze the 
renormalizability of the theory to all orders. Let us first identify the final complete action to start with. 
For this purpose, we need to introduce additional external sources 

{a"^,U' ,Y" Xi :^f) in order to 
define at quantum level the composite operators entering the nonlinear BRST transformations of the 
fields {A"^,L",B'f,B",Gf,G"), eqs ( |Z4l ). In the present case this term reads 

5ext = ^ J d\{-Q."^Al + l^"(:" -y'iB- - Y^^ +X"G'l +Xf'^i) , (2.36) 

with 

sQ-l = sL" = 0, (2.37) 
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X X Y Y D. L 


dimension 
ghost number 
Qf -charge 
2/' -charge 


3 3 3 3 3 4 
-2 -1 -1 -1 -2 

1-11-10 



Table 3: Quantum numbers of the external sources 



and 



sYf' 


— -^i ) 


sXf 


= 0, 




= -Y- 


sY'; 


= 0. 



(2.38) 

The quantum numbers of the external sources 

{Q."^,L",Yl,Yi",X''i,X^) are displayed in Table 12 

Furthermore, we have to add the extra source term Sextra for renormalization purposes, as it was 
explained in Il40ll45l 

yahcd ^ ^ , ^ 

H — -^{Bi^i —GiGi){BjBj — GjG j) + \BjG"Vif,yU jjj^ + GiG^jUi^U y^v 

+Bf B'jVi^yV jfjN — Gj B'jVj/jyUi^yj — G"B'jU i^yjV jfjy + G,- BjUi^^Vjij^ 
\ y Y ^ ^ \ ^ ^ 

— ^^"^"ViijyjVj^y + -G"G'^jUi,jyU j^M — -^^B^BjVi^MVj^y, + -Gf GjUi/jyUj^y 

-^{U ifjvUifjvU japU jap + Vi^vVifjyV yapVyctp - 2U i^vUifivV ja^Vjaf,)^ , (2.39) 

where , , Xi , X2 > ^ are free parameters, and the gauge invariant rank 4 tensor X,"^'^'^ has the following 
symmetry properties 

y^ahcd y^dab ybacd 



and it obeys a generalized Jacobi identity 

j^manymbcd _j_ jmbnyamcd jmcny^abtnd jmdny^abcn 

Thus, the complete action we are looking for is 

^ — ^1 ~l~ »^ext ~l~ •S'extra 



(2.40) 
(2.41) 



Sym + 1 A l^b" d.A^+rd.Dfc' + '^;!M - CO? 51^ "'cot + ^r''co?a^(^t<c^) 



+Xi {B"B" — Gj G")(y jfjyVjfjy — U jj_,\,U y^v) 
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I abed 
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+Bj ByV^^v^;>v — Gj B'jVjfj^^Uifjv — G"B"jU ,-^,vV y^v + Gj BjUi^vVjfj^ 

2 



-B"B"jVifjyV y^v + -jG"G"U ifjy,U jfjy, — -Bj BjVifjy/Vj/jy, + "G; G jUij_,yU y^y 



fubcv^ „b nc 



rabcvu „b'D'- 



rabcva „b 



(2.42) 



3 The complete set of Ward identities 

In this section, we have enlisted all known Ward identities, associated to the action (I2.42I ). 



The Slavnov-Taylor identity 



_6^6I 6181 81 ^81 -t, 8£ tt^ 8^ 
8il«8A« + 8L^8?^+ 5^^^''5^^'^'S^~ 
/ 8£ _8I 81 8£ 8£ 



81 



81 



^^8M^ 
81 



^ '^'8l^-,v 8X:' 8F,« 







(3.1) 



• The Landau gauge fixing 



5b" ~ 



(3.2) 



• The antighost equation 



8£ 8£ 

8^^ 



0. 



(3.3) 



The ghost equation 



5^ rabcf b^^ ./,5r .81 81 81 



8c' 



V 



class 



(3.4) 
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The rigid group transformations 



Vic = {A,b,c,c,(^,^,(0,(0,B,B,G,G,a,L, 
M,M,N,N,X,X,Y,Y) . 



The SL(2,M) invarianc^ 



The (t)-equation 



• The co-equation 



The (|)-equation 



The co-equation 



The global C/(/) invariance, / = 4(A^2 _ 



f ^4 /^^« 5£ 51 „ 51 ^ 51 51 



—a 51 51 —a 5£ 

2 See also (651. 
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• The rigid symmetry related to the horizon function 



• The symmetries relating the auxiliary fields (|), (|), co, co to the Faddeev-Popov ghost and antighost 

. /A(w§-c»§+r^^)^0. (3.>2, 



The global U{6) invariance 

SL 



Q^jiL) = J dx^B, — -Bj^ + G, — -Gj^ + U,^^-Uj^ 



SL — SL ^ SL —a SL ^ SL — « SL . 



(3.14) 



• The rigid symmetries related to the mass operator 

^^.J"^(E) = 0,aG {1,2,3,4} 
with 

(1). . r , / , 5L -aSL SL - 5L „ 6£ ^ 51 \ 

("3") /„N /" .4 /—a 5E — a 5E — SL — SL —a SL —a 5L \ 

(3.15) 



4 The linearized operators and (anti-)commutation relations 

In order to facilitate the upcoming vast amount of algebra required for the determination of the most 
general counterterm, we shall give her some (anti-)commutation relations between several (linearized) 
symmetry operators. 
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The equations (13.11) . (13.61) and (13.131) generate, respectively, the following linearized operators: 
61 6 6£ 5 5£ 6 61 5 ..5 „ 6 5 



r 4 52: 5 52: 5 52: 5 52: 5 5 „ 5 -^^ 5 

__ _5_ /_5r 5^ 8 



®Z = 



^4 , « 5 5£ 5 5£ 5 

T r 1 1 



8 5£ 5 _ 5 8 



8c" 8^^ 



8 8£ 8 81 8 



^^8^" 8L"8(0? 8(0^8L" 



0. 



Consequently, we are able to derive some useful (anti-)commutations relations: 

8 



8 

8^ ' 



_8_ 

8c"+ ^'8a^ 

0, 



_8_ 



_8_ _8 



. 5 













_8 
8^? 



5 ^ 

- -I- ;^ — o f'^^ 



'8F 



81 



8^^ ^dM'/j 



rate /^T^ SI _^ 8E 8E 



Q<7 
0, 



(4.1) 
(4.2) 
(4.3) 



{^.f' , .4 = /A (8.8, + 8.8,.) {Gt^-U,,.^^-X^^[ 



^8F/" 
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{^f , = -/A (5.5,-5.5,.) (^^J. -^Vvg|- + r,^' J.) 



5 



(4.4) 



5 Characterization of the most general counterterm 

In order to characterize the most general invariant counterterm which can be freely added to all or- 
ders of perturbation theory, we perturb the classical action £ by adding an arbitrary integrated local 
polynomial Zct in the fields and external sources of dimension bounded by four and with zero ghost 
number, and we require that the perturbed action (Z + riZcT) satisfies the same Ward identities as £ to 
the first order in the perturbation parameter r\. Making use of the BRST cohomological results L64il . 
we may write that 

rcT = ao5YM + ®2:A(-i) , (5.1) 



where S^; is the nilpotent Unearized Slavnov-Taylor operator of eq. (l4.1l ). 

SzSi: = . (5.2) 
The expression A^^'Ms an integrated polynomial of ghost number —1, in the present case given by 



A(-i) 



+017X3 ^GjG'jVi^yU jfjw + Gi^'jViijyV — -B'lB'jUi^vV + -B'lG'jUi/^ivU j^y 

Y Y ^ ^ \ y^abcd _j_ abed ^ , , 

- 2 V,,uv + -GiG j Vif_N t/i>(v j + a 1 8 G,- B,- [BjBj - GjGj j 

+ai9 X A^,M«/ + a2o gf'%j ^'lA^ + ^21 gf^M^j^A'^ + 022 gf^^A', d.^'j 
+a2,gf''Wi{d,Al)(Sf\ + a2,gf'''&'iA'^,d,B^^ 

+a2^gr''Ui,,A%B'i+a2igr'%,yA'^'yVl + ^^^^^ . (5.3) 

Due to the Ward identities given in section 4, the counterterm Zct must obey the following constraints 

SZSCT=0, !D5:rCT=0, ^"rcT=o, 'W^"rcT=0, 

(i-,EcT = o, ^.y'''''')rcT = o, K.f'\«^}rcT = o, ,%^''^}Ect = o, 

5ScT 5ZcT 5ZcT ^ 5ZcT 5Zct ^ 52Lct 5Zct ^ 
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5^? 
51 



8£ 



0, 



'SN' 



-0. 



By applying the constraints (I5.4I ). one can show that 



abed 



2 face _redb 



(5.4) 



(5.5) 



and 



+ [(ai + ^2 + «4)2a^A« + (2fli + a2 + a4)5r'''^A^^] • + ^^vB- + t/,>vG; - F,-^vB- ) 

{W^B^ - ^G^)(Byj - ^Gj) + (2a2 + a,)X, {W^Bt - ^G^) 



+{4a2 + a5[ 



yabcd _^ 

16~ 



: (V'/^vV^^jv — U ifjvUi/jy) + {laj + a-])Xj, yBj G^jVi^-^U y^v + G,- G"Ui^yU + B^ B"jVi^y,V y^v 



-G,- B'jVjij^Uifjv — G"B"jU i^^V jfjv + G,- B^Utfj^Vj^y, — -BfB"ViuyV /i,v + -G"G"U iuyll , 



--BlB^jVi^vVj^v + -Gi GjUi^vU j^y^ + as C,{U i^vUifj^U ja^U + V/^iv^/zv^ yap^yap 
-2?7 /^vt^i^v^yapV'yap) + Xl (^(>v9^V';^v - ifiv^^Ui^v) 

where we renamed the coefficients a„ as 



«3 

an 



«2, 

—2(34, 



an 



as, 
ay, 



ai2 
^13 

ai4 



«8, 
ag, 

aio. 



(5.6) 



(5.7) 



6 Stability of the action at the quantum level and renormalization fac- 
tors 



As a final step, we must show that the most general counterterm £ct can be reabsorbed by means of 
a multiplicative renormalization of the parameters, fields, and sources already present in the starting 
action £. Taking 



¥o 

^0 



1 abed 



^1/2 

ZjJ, 



2 yabed _j_ ^^abed 



(6.1) 
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where 



\\f = {A,b,c,c,<^,(^,03,03,B,B,G,G} , 

J = {^,L,M,M,N,N,U,TJ,V,V,X,X,Y,Y} , 

^ = {8,l,%i,%2,^,h,h} , (6.2) 



we must show that 



r(i|/o,7o,^o) = r(¥,-/,^)+iircT(¥,-/,^) + 0(ri2) • (6-3) 
After some algebra, for the renormalization factors {Z} we obtain 

Zh = , Zc = Zf = Zg z^ ^ , 

Z(|) = Z^ = Zg Z^ , Zq) = Zg , Z(j) = z^ , 

7 _ 7__7l/2_ 7-1/27-1/4 

7_ _ 7l/2_7-l 7 _7l/2_ 7-1/2 

7 _ 7I/2 _ 7-1/27-1/4 _7l/2 
— Zc — Zg Z^ ,Zt — Z^ , 

V V V '7^^/2^71/4^7-1/2 

^ _ — Zy — Zj7 — Zg Z^ Zg , 

Za = l+ri(ao + 2ai) , Zg = l-riy, 
= Zj = Zg = Zg^ = l+ri(2fl;2 + a3) , 
Zy = Z^r = Zf/ =Zj7= l-ri(^y + y-c?4 

Zx = l+Tl(a5-2a3), Z'^'^'' = ^{4a2 + a5H'''''-' , 

= l+r\iao-2a4 + a(,) , Zx^ = I + r\{aQ - 2a4 + aj) , 

Z^ = 1 +ri(2ao + 2a3 — 4a4 — ag) , Z,^, = 1 +ri(fl;o + «3 — 2a4 + a9) , 

Zj^j = i+r\{ao + ai,-2a4 + aio) , (6.4) 

hereby confirming the renormalizability to all orders of perturbation theory of the action ( I2.42I ). We 
draw attention to the fact that the renormalization of the quartic tensor coupling X"^"^^ involves an 
additional additive part given by z"'"^'^ 1.62] . As a consequence, X"^'^'^ = is a not a fixed point of 
the model. This originates from the fact that the interactions proportional to the other coupling 

reintroduce by quantum effects the tensor coupling oc [b"^^B^^^ — G^^yG'^y) ^Bpc^pa ~ ^pa^pcj ■ This 

is nicely reflected in the renormalization group function of X"^'^'^, that was calculated at one loop order 
in l62l using dimensional regularization ind = 4 — 2s dimensions and the MS scheme. 



OfJ 



12CAr''''^a + ^CAf^Pf^Pa^ + leCAf^Pf'^a^ + 96df"^a 



(6.5) 



where a = and d'^'^'^ is the totally symmetric rank four tensor defined by d'^'^'^ = Tr [T^Tj^T^Tf^ 
Clearly, we have 
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6.1 The physical action and some of its properties 



We have shown the renormaUzabiUty of the complete action ( I2.42I ). In particular, since the renormal- 
izability holds for all possible values of the sources, we have also proven it in the case that the external 
source part, ( |2.36b . is zero, while the other sources attain their physical values dictated by ( I2.17I) and 
(12.261) . yielding the complete physical action 



•^physical 



+ 

+ 
+ 
+ 
+ 



d\ 



-m 



1 abed , 
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d^x ( —^m^ 



(6.7) 



Apparently, the symmetry content of the action (12.421 ) given in section 3 is sufficiently powerful to 
avoid mixing between the Zwanziger fields/sources on one hand and the mass related fields/sources 
on the other hand. 

The term oc C^m^ in the final action (16.71 ) is irrelevant for the renormalization of Green functions, but it 
becomes important when one looks at the renormalization of the vacuum energy E{m). The parameter 
^ is the so-called LCO parameter, and its value ought to be fixed by requiring a homogenous linear 
renormalization group equation for E{m), whereby ^ is made a function of the available couplings. 
This point is however beyond the scope of this paper, the interested reader is kindly referred to e.g. 
BSl l46l l5n l66l for more details. It is a remarkable feature of the Zwanziger action that there is no 
need for such a LCO parameter in front of the y*-term in the action (16.71 ) [i40ll4n[T4ll . 

If we make abstract of the Gribov-Zwanziger part for the moment, we estabUshed a "supersymmetry" 
for the action S"^-^y^ll~^ , generated by IIH 



5. J 5! 



5? 



vii=0,y=0 
physical 



= 0, 
0, 



for all other fields »F , 

0, 

0. 



(6.8) 



We used this symmetry in 11621 to show that the massless version of our gauge model is equivalent 
with Yang-Mills ordinary gauge theories, despite the extra (quartic) interactions between the fields 



B^y, and G^y. A completely similar S^-cohomological argument as presented in [62,1 can 
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be used here to actually prove that the action ^physical original Gribov-Zwanziger action give 

rise to the same Green functions at any order of perturbation theory when we restricts ourselves to 
those functionals built from fields in the original Gribov-Zwanziger action, meaning that the quartic 
coupling X"''"' cancels out from the final results. 

The combination of the previous result and the already mentioned absence of mixing, also implies 
that the already known renormalization group functions and relations for the Gribov-Zwanziger action 
BOl [1411411 and massive gauge model fl31l62i remain valid when both are combined into one action, 
at least whenever massless renormalization schemes like the MS one are employed. 

When the sources are set equal to their physical values (12.171 ) and ( 12.261 ) in order to obtain the action 
S'physicah the BRST symmetry (12.41 ) is however broken. It is worth having a somewhat more detailed 
look at this. 



7 The breaking of the Slavnov-Taylor identity scrutinized 



7.1 The case of the massive gauge model without the Gribov restriction 

In order to avoid too lengthy expressions, we shall momentarily skip the Gribov restriction, and con- 
centrate on the massive gauge model already studied in earlier papers fl31l62l . 



Let £ thus be the complete action given by 



+ 



y abed 
1^ , 



(B-Bf - ^Gf)iByj - G^Gj) + X3iB"Gyi,,Uj,, + (^Gpi,^U 
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+Bj B'jV^ivV j/jy — Gi B'jVj/jvUi^v — G'j'B'jU i/jvVjjj^, + Gj ByL^i^v Vy^v — B'jV i/j^V j/jv 

i2,, 77 :i2, 



This action E obeys the Slavnov-Taylor identity 
with 

51 '~ 



(7.1) 



(7.2) 



_5^5£ 5£5£ 51 /_5£ p« ] ^ 5151 5£ 5£ 

U^h^^UJ'^'^ 5i^^ I 5Ff ^ j 5Bf ^ WfWi ^ SX" 5Gf 



+ 



51 5£ - 



51 5i: 



I/JV 



—aSL 5£ 
'5r 5F,« 



(7.3) 
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Since the theory is stable and free from anomalies at the quantum level, we may write down a renor- 
malized IPI quantum vertex functional f64], 

r = i+nr (7.4) 

which fulfills the quantum version of the Slavnov-Taylor identity ( I7.2I ). i.e. 

5(r)=0, (7.5) 



Sir) 



d X 



6r 6r 6r 6r 



5r 



6r 



6r 5r 6r 5r 5r 



6c' 



+ 



5r -a\ 5r _ 5r 



5r 



6B« 57/' 5B; 6x;'5G: 



8X/' 



6G' 



517 



8Vi 



8r 



(7.6) 



Let us now analyse the quantum properties of the model when the sources attain their physical values 
( 12.261 ). First of all, let us give a look at the classical action Zp/,, obtained by £ by setting the sources 
to their physical values, namely 



physical value of (Vjfiy .V jf^, ,Ui^ni ,U 



or explicitly 

£ph 



+ / d^x (b"d,A" + c"d,Dfc'' 



•\ab r\bc nc 



+ 



d'^X 



d'x 



rfl Tib f^" f^b 



bLbL - g'g'! 



pa" pa 



'pa^pa 



d X I —Lm 
1 4^ 



(7.7) 



(7.8) 



It is easy to check that Zph is not BRST invariant w.r.t. (12.41) . In fact, it turns out that 

,2 



im 

^Sph = — 



m 
16 



d^xlB 



(7.9) 



This equation shows that the breaking of the BRST symmetry (12.41) is not linear in the quantum fields, 
and hence the breaking terms have to be treated as composite operators [64|. Therefore, equation (17.91 ) 
cannot be renormalized as it stands. The two breaking terms have to be taken into proper account. 
This is precisely achieved by introducing the local sources (V,y,v5 V',y,v,?/i^v;f^!>iv)- In other words, 
these sources allow us to take into account e.g. the presence of / d^xG^^F^^ and its renormalization, 
which is expressed by the renormalization factor of the source Ui^^j. 

We would like to understand what happens to the BRST symmetry at the quantum level, when the 
sources attain their physical values. It is instructive to study this limit by means of the Slavnov-Taylor 
identity (17.51 ). Let Fph be the IPI functional obtained from F when the sources {Vt^^ ,V i/jv ,Uifjv ,U i/jv) 
attain their physical values 



ph 



physical value of {Vjf^ ,V j^i^ ,Uif,v ,U if^) 



(7.10) 
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We can write 



d xViijy-T= — 



im 



physical value 



m 



physical value 



(7.11) 



- physical value 

Stands for the generator of the IPI Green functions with the insertion 



where e.g. [{J d^xG^^^F"^) -F 

of the composite operator (/ d^xG"^F") . Of course, it holds that . . . • F 

L J physical value 

follows that the quantum action Fph obeys the broken Slavnov-Taylor identity 



..r 



ph 



It 



'^(rph) = ^ 



d^xG%F^^-r 



m 



physical value 



d^x [B"^y - B^v) • r 



, (7.12) 



physical value 



where 



^(rph) 



d X 



.5r, 



ph 



5a« 5A« 



5L" he" 



8c" 



+ 



6r, 



ph 



5F" 



+ G" 



SFph ^ SFph SFph 



5Bf 57." 



+ 



sr, 



ph 



5Xf 



sr, 



- ph T?-u pii 

=7) ; ~r A; 



«5r, 



ph 



,5r, 



ph 



5G" 



6xr 



6F" 



(7.13) 



It is worth underlining here that the equation (17.121 ) is in fact nothing more than a direct consequence 
of the Slavnov-Taylor identity (17. 5I ). when the local sources (V/^y > ^i>jv , ^^r>jv > ^^i>jv) attain their physical 
values (12.261 ). 



We conclude that Fph does not obey an exact Slavnov-Taylor identity. Of course, (17.121 ) translates at 
the quantum level the fact that the classical action Zph, obtained from £ by bringing the sources to 
their physical values, is not BRST invariant, according to (17.91 ). However, even if Fph does not obey an 
exact Slavnov-Taylor identity, (17.121 ) has far reaching consequences on the behavior of the IPI Green 
functions obtained from Fph, i.e. when the sources are set to their physical values. 

Let us consider the breaking term [ ( / d'^xG^^F^^) ■ F] . Typically, Slavnov-Taylor identities at the level 

of Green functions are obtained by acting with a test operator like S(p(a-i) S(p(a- ) ' "^^^^ ^ generic 
field, on expression (17.121) . and setting all sources and fields equal to zero at the end. The condition to 
be fulfilled so that the breaking would be harmless is quite easily found, since the breaking term will 
vanish whenever 



5"[(/AGgyF;y)-F] 
8(|)(xi)...8(|)(jc„) 



AG^y(x)C(^) )(p(xi)...(p(x„ 



(7.14) 



IPI 



meaning that the IPI Green function with the insertion of the operator (/ d^xG'^^F^^) and with n 
amputated external (|)-legs should vanish. Thus, if the condition (17.141 ) holds and an analogous one 
for the other breaking term, the right hand side of (17.121 ) is harmless, so that everything goes as if the 
theory would fulfill an unbroken Slavnov-Taylor identity, namely 



5(Fph)=0. 



(7.15) 



24 



The set of identities for which this happens is quite large. Certainly, it contains all Slavnov-Taylor 
identities which are obtained from (17.121 ) by acting only on the original Yang-Mills fields or even the 
B"^^ and B^^ fields. In this case, there is no way to obtain a nonvanishing contribution to the breaking 
term because of the presence of the G^^-ghost field in the right hand side of (17.141 ). For example, the 
Slavnov-Taylor identity for the IPI gluon propagator can be obtained from ( 17.121) by acting on it with 
the test operator Sc(,.)sa (y) ^^'^ setting all fields and other external sources equal to zero. The breaking 
terms will be irrelevant as the Green function ((/ t/^zG^y(z)F^" (z)) c(x)Av(3')) as well as the other 
one are trivially zero. 

We conclude that most Green functions will behave as if the theory obeys the unbroken Slavnov-Taylor 
identity ( 17.151 ). The same considerations outlined for the Slavnov-Taylor identity can be repeated for 
the other Ward identities. The corresponding breaking terms will always contain the integrated ghost 
fields G|[5v and/or G^^ which, in most cases, will lead to vanishing contributions when inserted into 
a Green function, thereby making the breaking harmless. The beauty in all this is exactly the fact 
that the breaking of the Slavnov-Taylor and other Ward identities can be brought under control at the 
quantum level by the introduction of a suitable set of local sources. All the renormalization results 
of the action with arbitrary values of the sources are then preserved once the sources are put equal to 
specific values. 

Since the classical part of the action (17.81) . obtained by skipping the gauge fixing term (11.101 ). is gauge 
invariant w.r.t. the gauge transformations ( 12.241 ). we expect that there should be a nilpotent BRST 
generator at the quantum level for the gauge fixed action (17.81 ). Nevertheless, we have just seen that 
the BRST operator (12.41) no longer generates an exact symmetry of the action (17.81 ). As it was already 
discussed in 11621 . the following nilpotent transformation 





= -Dfc\ 


s'c" 


S rate a b 
- c C , 


s'B" 

"/JV 


= gf'c'B'^, 


sB^^ 


= gf''c%. 


s'G" 


= gf'c'Gl, 




= gf''c'G% 




= b", 


s'b" 


= 0, 


s 


= 0. 



(7.16) 

generates an invariance of (17.81) . One shall easily recognize that there is an intimate connection be- 
tween the transformations s (12.4b . s' (17.161 ) and 5^ (16.81) . namely we have 

s = s' + 5,. (7.17) 

Then we can say that the breaking of the BRST symmetry s and its associated Slavnov-Taylor identity 
is entirely due to the loss of the supersymmetry 8^ when the physical limit ( 12.261) is taken [62|. 

Evidently, the unbroken BRST symmetry s' can be used to construct unbroken Slavnov-Taylor identi- 
ties between the Green functions of the massive gauge model ( 17.81 ). 
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7.2 The case of the massive gauge model with Gribov restriction 

A very similar analysis can be made when we consider the full action ( I6.7I ). In this case, there are ad- 
ditional breaking terms coming from the physical limit (12.17b of the Zwanziger sources. In particular, 
applying the transformations s or s' on (16.71 ). we find 



Irrespective of the choice of BRST transformation s or s', the physical Gribov-Zwanziger action is not 
BRST invariant anymore. However, repeating the same argument given in the previous subsection, 
it turns out that the breaking terms can be treated consistently at the quantum level, leading to a 
renormalized broken Slavnov-Taylor identity. Furthermore, in most cases, everything will go as if the 
theory obeys an unbroken Slavnov-Taylor identity, because the breaking terms are in fact harmless. 

8 Conclusions 

In this paper, we have shown that the nonlocal gauge invariant operator Tr / d^xFfj^,{D^)^^ F^j^ can be 
coupled to the Gribov-Zwanziger action in a localized form. By embedding this model into a larger 
class of models with local sources, we established a comprehensive set of Ward identities, which were 
sufficient to prove the renormalizability to all orders of perturbation theory. Specializing thus to a 
particular values of the local sources, we conclude that we have constructed a renormalizable action 
( 16.71 ) that allows us to study the gauge invariant mass m in combination with the restriction to the first 
Gribov horizon, obtained when the effective action is minimized w.r.t. the Gribov mass parameter j^. 
This restriction gives a first source of nonperturbative effects in gauge theories, as explained in the 
introduction: the gluon/ghost propagator gets infrared suppressed/enhanced, while the Gribov mass 
y is fixed in terms of the QCD scale Aqcd by means of the requirement that the effective action is 
minimized with respect to it. 

We also payed attention to the breaking of the BRST invariance when the sources are set equal to 
their physical values. We have elaborated on the fact that in most cases, the breaking terms in the 
Slavnov-Taylor identity are harmless, since they usually induces zero contributions to the identities 
between Green functions. 

In the main body of this paper, we have extensively used and studied the BRST invariance in relation to 
the renormalizability. However, there is another major reason why the BRST symmetry is so important 
for perturbatively handled gauge theories. In a certain sense, the BRST invariance is the quantum 
version of the gauge invariance, and as such it should play a major role in reducing the number of 
relevant (physical) degrees of freedom, at least at the perturbative level. It is well known that a BRST 
symmetry with corresponding nilpotent charge is a very powerful tool in establishing the unitarity of 
gauge theories at the quantum level once a gauge has been chosen, see e.g. |[37ll67ll68l 16911701 . 

When we discard the Gribov restriction, we retrieve the gauge model studied in fl31 16211701 . enjoying 
the BRST symmetry with nilpotent generator (17.161 ). Therefore, hope existed that the model might be 




(7.18) 



or 
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unitary, i.e. that one would be able to define a physical subspace :?/phys of the total Hilbert state space 
.?/ , such that ^pbys is endowed with a positive norm. This optimism turned out to be flawed, as it was 
shown in [|70J that the massive gauge model is not unitary. 

In addition to this, the restriction to the Gribov horizon only makes things worse. First of all, we have 
lost the (nilpotent) BRST symmetry, so any potential discussion of the unitarity cannot be based on 
BRST related tools. Further, we already mentioned in the introduction that the Gribov restriction gives 
rise to an infrared suppressed gluon propagator, and this suppression is so that the gluon propagator 
shows a violation of spectral positivity. Hence, the gluon is not expected to represent a physical 
particle, implying that we should certainly not expect unitarity from the Gribov-Zwanziger action 
when the gluons are treated as physical particles. We should rather expect the opposite. A hint 
that the Gribov restriction destabilizes the gluon is also given when we take a look at the tree level 
propagator, which in our conventions is given by 1511 



with the gluon form factor 



which can also be written in a "standard" propagator form 



'DiP^)= 4 , o ?.r.A > (8-1) 



'D(p^) = - +- , (8.2) 

i.e. as the sum of 2 propagators with imaginary masses squared. 

A similar reasoning can be applied when we do not implement the Gribov restriction. As we already 
outlined in 1 70], we can see the massless version of our model (17.81) as an alternative to ordinary Yang- 
Mills theory at high energies, based on their equivalence in the perturbative region [62]. The benefit 
of using our gauge model is that it is possible to couple mass terms oc m to it without jeopardizing 
the renormalizability. Then one can start looking for a sensible gap equation in order to generate a 
nonperturbative mass scale m. Said otherwise, we could start looking for a gauge invariant dynamical 
mass generation mechanism. The generation of such a mass parameter would break the unitarity at 
the level of the gluons, but we must recall that unitarity is only a prerequisite for the physical degrees 
of freedom. We can depart our research from the massless (unitary) theory |70]| . but we are no longer 
interested in describing the perturbative asymptotic high energy regime of QCD, but instead we are 
entering a phenomenologically interesting region where e.g. the gluons already loose their physical 
meaning as observables. In this energy region, the gluons should be rather seen as a kind of quasi 
particles with a finite lifetime, thus not entering the asymptotic physical spectrum, which we do not 
know how to describe. We can continue to use the gluon propagator etc., albeit the versions corrected 
by nonperturbative effects, like an effective gluon mass and/or Gribov restriction. At first instance, 
we can concentrate on the case with only the mass m to be fixed, but at a later stage, we can also study 
the influence of the restriction to the Gribov region Q., since we have just proven the renormalizability 
of this action to all orders. 
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